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ABSTRACT 

Lyct line emission from the Loeb-Rybicki (LR) halo, which is the expanding 
HI IGM (intergalactic medium) around the first star clusters and the ionized in- 
terstellar medium, is investigated by solving a radiative transfer problem. While 
the initial scattering optical depth is ~ 10 5 for the Lya photons, most of the Lya 
photons can escape when the cumulative frequency-shift due to the expansion of 
the HI IGM becomes significantly large. The current paper improves upon previ- 
ous treatments of the scattering processes and the opacity for the Lya transfer. 
Confirming the previous results of the LR halo, we investigate the effect of the 
aspherical expansion of the IGM. The asphericity is hypothesized to follow the 
initial stage of the gravitational deceleration to form the large scale filamentary 
structure of the Universe. According to our results, the effect of the aspheric- 
ity lets the peak wavelength of the line profile shift to longer wavelengths and 
the FWHM of the profile become wider than those of the spherically expanding 
model. To detect these features is meaningful if we are interested in the initial 
evolution of the large scale structure, since they reflect the dynamical properties 
of the IGM at that time. Furthermore, given the recent discovery of the high red- 
shift cosmological reionization, we briefly comment on the effects of the redshift 
and the cosmological parameters on the line profile. 

Subject headings: line: profiles — radiative transfer — galaxies: intergalactic 
medium — galaxies: formation 

1. INTRODUCTION 

The epoch of cosmic reionization is one of the most important epochs to understand in 
the investigation of the history of the universe, because it contains the imprint of important 
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information on cosmic structure formation and evolution. This is because the reionization 
of the intergalactic medium (IGM) occurs promptly (e.g. Gnedin & Ostriker 1997) and the 
cosmic ionizing background radiation can penetrate the dense HI condensations of virialized 
objects (Miralda-Escude, Haehnelt, & Rees 2000; Barkana & Loeb 1999). 

The epoch of the reionization observationally relates to the Gunn-Peterson trough (Gunn 
& Peterson 1965). The significantly ionized IGM in the observed spectra of QSOs at z < 5.6 
(Dey et al. 1998; Hu, Cowie, & McMahon 1998; Spinrad et al. 1998; Weymann et al. 1998) 
implies that most of the intergalactic hydrogen in the Universe has been reionized before 
z ~ 5.6. Some standard cosmological models predict that reionization takes place around 
a redshift of z ~ 6. However, according to the results of the WMAP project, Kogut et 
al. (2003) may have discovered the high redshift reionization of the Universe (z ~ 20). To 
resolve this discrepancy, the direct detection of the emission from the first objects is very 
important. This is because that emission reflects the physical and chemical properties of 
the first star clusters (Shchekinov 1991; Kamaya & Silk 2002) and/or protogalaxies. Once 
we know the origin of the first emission, the ionization rate at that time is investigated. 
Therefore, new detailed studies and confirmation of the previous studies of the first emission 
in another way are important for the subsequent progress of the research fields. 

It is interesting to detect just the first Lya photons since the star formation rate and/or 
the ionization structure of the interstellar medium are examined from that intensity. Using 
standard cosmological models, cosmic recombination takes place at z ~ 10 3 . After that 
epoch, the Universe becomes predominantly neutral. When the first galaxies and/or star 
clusters form, they are surrounded by the neutral IGM. One might naively expect the damp- 
ing wing of the Lja trough to eliminate any trace of a Lyct emission line in the observed 
source spectra. However, the Lya photons are scattered and diffused redward in frequency 
owing to the Hubble expansion of the surrounding HI IGM (Loeb & Rybicki 1999 ; hereafter 
LR99). When the cumulative frequency- shift becomes significantly large, they can escape 
and travel freely toward the observer. 

Although we investigate a similar problem to that in LR99, the radiative transfer model 
is improved in the two ways : (1) Formulation of the geometrical scattering (i.e. test photons 
scatter spatially three dimension); and (2) Formulation of the Lya opacity. Thanks to these 
two improvements and the confirmation of the previous studies, we can study more complex 
physical conditions in the expanding HI gas. Since the large scale structure (LSS) forms 
filaments, this suggests the initial stage of the gravitational deceleration for the LSS favors 
an asymmetrical expansion of the HI halo. Indeed, the inhomogeneity of the self-gravitating 
system can grow in a free-fall time (e.g. Lin, Mestel, & Shu 1965). As long as the expansion 
time-scale is comparable to the free-fall time scale, the asymmetrical deceleration of the 
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IGM occurs in the time-scale of the Hubble expansion, although the initial asymmetry is 
obviously very weak. As the first and the simplest hypothesis to investigate the effect of 
the inhomogeneity via the gravitational deceleration for the LSS formation, we adopt the 
aspherical expansion model of HI IGM. The physical implication drawn from the current 
analysis will suggest to us further research projects for the first Lya. Since the lumpiness 
and/or small-scale irregularity may be important for a realistic treatment of the first Lya, 
we will investigate these effects in subsequent papers. 

This paper is organized as follows: In §2 we describe the formalism of the Lya line trans- 
fer problem commenting on the differences between our approach and LR99. In §3 we outline 
our Monte Carlo simulation method. In §4 we discuss our numerical results including the 
effects of a more complete representation of the Lya opacity, the effect of asphericity in the 
expanding IGM, and the sensitivity of our results to redshift and cosmological parameters. 
In §5 we summarize our study. 



We formulate the problem according to LR99, while the HI IGM is hypothesized to 
expand aspherically. In general, the exactly spherical expansion is not expected as long 
as the IGM decelerates to form LSS. This is recognized because the free-fall time, is 
comparable to the Hubble expansion time, th, as explained in §1. Then, to find the effect of 
the asymmetry, we shall set the most simple assumption to the recession velocity, spheroidal 
expansion: v = H s [e 2 (x 2 + y 2 ) + z 2 ] 1 ^ 2 , where v is recession velocity, H s is the Hubble 
expansion rate at the source redshift of z s , and e is ellipticity or decreasing rate of the 
recession velocity to the pure Hubble expansion. We shall estimate the allowed order of e. 
The ratio of the two time-scales is denoted as 



Here, 5 is the ratio of the density of the halo to the mean density of the Universe at z. Of 
course, this ratio must not be different from unity very much as long as IGM expands. Since 
the density of the aspherically decelerating HI halo is larger than that of the spherically 
expanding halo, we find 5 — 1/e with mass conservation. Fixing the ratio of Eq.(l) to be 
unity and postulating 5 ~ 6, for example, we estimate e ~ 0.2 for VLm = 0.4. Thus, we 
recognize that order of e is about 0.1. The case of e — 0.1 is considered as an extreme case. 
Obviously, e — 1 corresponds to LR99. 



2. 



FORMULATION 




(1) 



Although e should approach to unity differentially from small scale to large scale, we 
assign e a single value for simplicity. Fortunately, the size of the concerned region, ~ r*, 
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which is a characteristic distance from the Lya photon source as defined later, is ~ Mpc as 
estimated in LR99. That is, r* corresponds to the size of groups of galaxies and the region at 
can be regarded as decelerating region although the allowed deceleration should be small 
(but not zero). Expansion law of the spheroid approaches to that of the spherical model at 
larger scale than r*. Properties of test photons are also determined at ~ r*. Hence, the single 
value hypothesis for e can present the physical insight into the effect of the deceleration for 
the LR halo. In realistic condition, by the way, e inside r* should decrease toward the center 
differentially. Then, our finding characteristic effect of the deceleration may be enhanced 
because the number of scattering can increase owing to the deceleration as shown in this 
paper. 

Let I = I(v,Q,,r) be the specific intensity (in photons cm -2 s _1 sr _1 Hz -1 ), where v 
is comoving frequency, r = (x, y, z) is the radial vector from the source galaxy, and Q, is 
the scattering direction. The spatial components of the direction cosines of f2 are (u,v,w), 
respectively. With the assumption of isotropic coherent scattering, the comoving transfer 
equation for a resonance line is given by 

dl 

n-VI + at/j(e,w) — = X v(J-I) + S. (2) 

Here, v is the shift of the frequency u a — z/ p hoton where v a is the resonant frequency and z/ p hoton 
is the photon frequency; \v is the scattering opacity at v (eq. (4) in LR99); J = (1/2) Id\i 
is the mean intensity; S is the emission function (in photons cm -3 s _1 sr" 1 Hz -1 , eq. (3) 
in LR99); a = H s v a /c is the increasing rate of the frequency shift per unit distance in the 
spherically expanding IGM (eq. (2) in LR99); and 

1>(e,w)= [e 2 + {1 - e 2 ) w 2 } 1/2 (3) 
represents the effect of aspherical expansion. 

In LR99 the scattering opacity \v wa s approximated as \v ~ fl/v 2 - In the current 
paper, we formulate the opacity in the following more complete form: 

X, - a^^<< W , (4) 

for v 2 ~ A 2 J16tt 2 , (5) 



u 2 + A 2 J16n 2 



^ for A 2 /16tt 2 < v 2 & v < u a , (6) 



1--J ^ for v~v a , (7) 
for v > v a , (8) 
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where v a = 2.47 x 10 15 Hz is the Lya frequency, A a = (8n 2 e 2 f a /3m e c\ 2 a ) = 6.26 x 10 s s" 1 
is the rate of spontaneous radiative decay from 2p to Is energy levels of hydrogen, f a = 
0.4162 is the oscillator strength, \ a = 1216 A is the wavelength of the Lya line, and (3 
is 1.5fib^o(l + -2s) 3 cm -1 Hz 2 (eq. (5) in LR99). In the initial scattering phase, Eq.(4) is 
important since the frequency shift can be smaller than A a /47r. When the photon diffuses 
and the frequency shift is comparable to v a , Eq.(7) must be adopted as \v 

We normalize the frequency shift and radius by convenient quantities like LR99. An 
appropriate frequency scale, v*, is estimated from t„ ~ 1. Here, r u is the Lya optical depth 
from the source to the observer, and is defined as 

POD 

Tu = / Xv+aT]>rdr. (9) 

Jo 

The characteristic frequency shift of is measured along the direction of z (i.e. w — 1). 
Eq. (9) with Eq.(7) can be described as 

rio / _ v _ 1± ar\ 4 g/g = 
Jo V ^ J (^ + mf ' 

where Iq = (f Q — is*) /a is a characteristic distance for the photon frequency to reach zero. We 
note that always depends on (3 and e. The appropriate radius, where the frequency-shift 
due to the Hubble expansion produces the frequency shift of v*, is given by r* = v^ja. For 
example, with e — 1, we find v* = 9.64 x 10 12 Hz and r* = 0.77 Mpc at z s = 10 with a 
set of cosmological parameters of Q b = 0.05, M = 0.4, D,\ = 0.6, and h = 0.65. 1 The 
normalized frequency shift v and distance f are then defined to be v = z//z/* and f = r/r*. 

For the reader's convenience, we present the Lymana opacity cm" 1 in Fig.l as a 
function of v. The solid line denotes our Xu, and the dashed line denotes the previously 
adopted one from LR99. According to the definition of z>, the photon frequency z/ photon 
becomes small as v increases and reaches zero at v = v a jv* ~ 2.56 x 10 2 . At v = 0, the 
frequency of the photon is the Lya frequency u a . The differences depicted in the figure 1 
are reflected in our results, while they are usually small. 



1 We adopt these cosmological parameters because of the comparison to LR99. Although the recent results 
of WMAP present slightly different cosmological parameters (Spergel et al. 2003), fortunately our results 
are not altered very much as shown in §4.4. 
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Fig. 1. — Lja opacities are presented. Solid line represents the fully formulated opacity 
[Eq.(4) in LR99], and dashed line does the approximated opacity in LR99. 
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3. Algorithm of the Calculation 

LR99 solved a dimensionless transfer equation in the diffusion regime, which is valid in 
the region v ^> f. They also studied the numerical solution of the transfer equation using 
a Monte Carlo technique and found that it agreed with the diffusion equation in the low 
frequency shift region. We also solve the transfer equation numerically by using a Monte 
Carlo method very similar to LR99. The number of test photons is 10 8 , the same as in LR99, 
and permits us to compare our results directly to theirs. During the calculation two photon 
decay is always neglected since it is difficult to destroy the Lyman a photons transferred in 
the surrounding HI. This is because the event rate of two photon decay is determined by the 
process of a transition between 2s - 2p owing to the radiative collision with an electron (Breit 
& Teller 1940; Spitzer & Greenstein 1951). As long as the relic electron number density is 
very small, the event rate of two photon decay should be negligible. 

LR99 adopted the initial shift method in the first step of each simulation. That is, the 
initial shift of v is 0.01 and f is determined stochastically by means of the diffusion solution 
of LR99. We have confirmed independently that the initial shift method is very reasonable. 
This is because the initial optical depth is always much larger than unity. Hence we also 
adopt that method in all first time steps of the simulations. 

1. 3D scattering 

LR99 assumed isotropic, coherent scattering through a spherically symmetric, neutral 
medium such that the scattering could be described by a single angular parameter 
fi = cos 6 relative to the radius r. We generalize to three dimensional scattering. The 
scattering processes in 3D coordinates is described by two parameters of fi = 1R\ — 1 
and = 7r (2i?2 — 1)- Here, is the change in azimuth, and R\ and R2 are random 
numbers on the interval of (0, 1). Defining fi as cos a, we measure a in the plane 
determined by the old and new photon paths. 

The direction cosines of u i+ i,Vi + \,Wi + i after the scattering are determined by using 
the direction cosines of the incident photon, Ui,Vi,Wi\ 

u>i+i = sin a (cos 4>WiUi — sin 4>Vi) / (l — wf) ^ 2 + cos aui, 

1 /9 

v i+ i = sina (cos^Wjfj + sin0Mj) / (l — wf) + cosaVi, (11) 
w i+ i = — sin a cos (l — wf) 1 ^ 2 + cosa-uv 

As \wi\ approaches unity, according to Witt (1977), the set of equations (11) breaks 
down and is replaced by 

Ui + i = sin a cos 0, Vi+\ = sin a sin 0, Wi + i = cos awi for 1 — |u>j| < 10~ 10 . (12) 
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2. Opacity selection 

LR99 used a single approximate form for the opacity as a function of frequency (i.e. our 
equation (6) for e — 1). We use the more complete form for the opacity by Equations 
(4-8). Selecting the appropriate opacity in Eqs.(4-8) according to v, we evaluate the 
path length along the ray, /, and the maximum optical depth: r max = J °° Xu+aipidl. 
Once the optical depth r for a test photon becomes larger than the maximum optical 
depth r max , we say that it escapes from the LR halo. Here, r is determined as — lni? sca 
where i? sca is a random number between zero and one. 

3. Aspherical expansion 

LR99 assumed that the neutral IGM surrounding the Lya point source undergoes uni- 
form Hubble expansion. Locally we do not expect this to be the case. The cumulative 
frequency-shift and r max differ depending on the scattering direction. According to 
Eq.(3), the increment of the frequency shift of the scattering photon is represented by 

= i>i + ip (e, Wi) I. (13) 
4. Results 

We first compare the results of our 3D simulation technique and improved approximation 
to the Lya opacity to those of LR99. We then expand upon that previous work by considering 
the aspherical expansion of the IGM. Recently Spergel et al. (2003) have used WMAP data 
to determine more precise values for the cosmological model parameters. WMAP data also 
suggest that the epoch of reionization may be at a redshift larger than 10 (Kogut et al. 
2003). Thus we conclude with a discussion of the effect of redshift and the cosmological 
model parameters on the observable Lya line profile. 

All figures presenting the results are made by counting the number of escape photons 
in the same interval of normalized frequency shift of v in log scale. We choose the interval 
of 10~ 3 from logz> = —1.0 to 2.0. For a clear presentation of the effect of the aspherical 
expansion, we count the photons all over the complete solid angle of 47r. 

4.1. Geometrical effect 

In Fig. 2 the geometrical effect of using our 3D numerical method for the scattering is 
compared to the 2D treatment of LR99. The same scattering coefficient used in LR99 is 
adopted with e — 1. The horizontal axis is the normalized frequency, i>, and the vertical axis 
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is the escape fraction of test photons. The 2D treatment of scattering by LR99 is the solid 
line; while our 3D treatment of the scattering is the dashed line. We observe no difference 
between the lines. This agreement comes from a simple approximation that the neutral IGM 
around the Lja source has spherical symmetry. Furthermore, it shows our 3D scheme has 
sufficient resolution to compare directly with the results of LR99. Hereafter, we compare 
our results with the spherically symmetric case in 3D coordinates, which is called 3D LR99. 



1 .0e-03 




100 



Fig. 2. — Confirmation of our scheme. Solid line is the spatially 2D scattering model and 
dashed line is the spatially 3D scattering model. 



4.2. Opacity evaluation effect 

How about the effect of opacity? Fig. 3 shows the results with the fully formulated 
opacity and e — 1 (solid line) and that of 3D LR99 (dashed line), using the same definition 
of both coordinates as in figure 2. It is found that there are almost no difference between the 
line profiles of ours and LR99's. This result confirms that the Lja opacity is approximated 
well in LR99. 




Fig. 3. — The effect of the selection of opacity is depicted. Fully formulated opacity is 
presented as the solid line (e = 1) and 3D LR99 is the dashed line. 



- 11 - 

Transforming v to the observational wavelength of A Q bs by 

(14) 

we compare our model with 3D LR99. In Eq.(14), the term [1 — (z^/V«) v] -1 represents 
the effect of the cumulative frequency- shift. The bin of the wavelength for photon counts 
corresponds to that of the frequency. We show the normalized luminosity of the Lya line 
per unit wavelength, L a (X), in Fig.4. This can be converted into the observational L a (X) 
by multiplying the steady emission rate of Lya photons at the source, N a (in photons s _1 ), 
by the Lja photon energy, hu a . The observed Lya wavelength emitted at z s = 10 is A Q bs = 

1 .2e-04 
1 .0e-04 
8.0e-05 

"i 

^ 6.0e-05 
4.0e-05 
2.0e-05 



1.33 1.34 1.35 1.36 1.37 1.38 1.39 1.4 

\ obs [ju m] 

Fig. 4. — Normalized luminosities of the Lya line profile as a function of observational 
wavelength A Q b s [/-tm] are depicted. Our result is solid line (e = 1) and 3D LR99 is dashed 
line. 

1.337 /j,m. As shown in the figure, the peaks of the line profiles shift to A pca k ~ 1.339 /xm. 
This shift from the line center is ~ 2 nm. We need a wavelength resolution of R ~ 7.8 x 10 2 
if it is to be detected. The full widths at half maximum (FWHM) of our line profile and 3D 
LR99 are 5.27 x 10" 3 /iin and 4.94 x 1(T 3 /mi, respectively 



Aobs — (1 + Z s ) \ a 



1 - 


^1 


1* 
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4.3. Aspherical expansion effect 

In this section the effect of the aspherical expansion of the HI IGM is summarized. Fig. 5 
presents the normalized luminosities of Lya line profiles for various values of the parameter 
e that characterizes the asphericity of the IGM. Here, the solid line is e = 1.0, long dashed 
line is e = 0.9, short dashed line is e = 0.7, dotted line is e = 0.5, long dash-dotted line 
is e = 0.3, and short dash-dotted line is e = 0.1. What is the distinct property of the line 

1 .2e-04 
1 .0e-04 
8.0e-05 
^ 6.0e-05 
4.0e-05 
2.0e-05 



1.34 1.36 1.38 1.40 1.42 1.44 

kb., [jii m] 

Fig. 5. — The effect of the aspherical expansion is presented. Each of the lines shows e = 1.0 
(solid line; e = 0.9 (long dashed line), e = 0.7 (short dashed line), e = 0.5 (dotted line), 
e = 0.3 (long dash-dotted line), and e = 0.1 (short dash-dotted line), respectively. 

profile with e ^ 1.0? The characteristic quantities are summarized in Table 1, where some 
results of other e's are also presented. The last column in Table 1 is the necessary resolution 
to distinguish A pea k from the Lja line center of (1 + z s )\ a . 

According to table 1, the distinct effects of the asphericity are the shift of the peak 
and the broadening of the FWHM of the line profiles. These are clearly seen in the line 
profile with e — 0.1. If there is strong asphericity of the expanding HI IGM, we will detect 
broadening of FWHM and the shift of the peak wavelength of the line profile. This could be 
observable with the JWST and signal the presence of the LR halo. 
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Table 1: Summary of the effect of the aspherical deceleration. 

For increasing asphericity, i.e. a decreasing e, the optical depth for photons to scatter 
becomes larger than in the spherical model. Thus, although the typical shift per scattering 
due to the expansion of the universe becomes small, the number of scattering events increases, 
thus increasing the probability that a photon will be shifted far from the line center. 

4.4. The Effect of Redshift and Cosmological Parameters 

WMAP data suggest that cosmic reionization may take place at z s ~ 20 (Kogut et al. 
2003). In Fig.6 we show the results of Lya line transfer on the Lya line profiles from a source 
at z s = 20 (solid line) compared to a source at z s = 10 (dashed line). The coordinates are 
the same as in Fig.3. We find the line profile of z s = 20 coincides with that of z s = 10 very 
well. This comes from the fact that most of the Lya photons can escape from the IGM when 
their cumulative frequency shift becomes ~ v*, and it takes place at the distance of ~ r* 
from the source. 

While V* and r* depend upon redshift, the shape of the two line profiles as functions 
of the normalized v are very similar. The effect of redshift is not significant except in the 
normalizations. This result presents the useful implication that when we examine line profiles 
at a variety of redshifts, it is sufficient for us to know the characteristic frequency shift and 
radius at each redshift. Hence, we summarize the redshift dependence of in Fig. 7. Here, 
the solid line is of Eq.(10), while the dashed line is z/^ pp which is adopted in LR99 and 
determined by Eq.(6) and Eq.(9). For simplicity, we present the cases with e — 1. We hope 
these results are useful for the heterorelation of future observations. 
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100 



V 



Fig. 6. — Line profiles of z s = 20 (solid line) and z s = 10 (dashed line) with e — 1 are 
presented. 




Fig. 7. — Redshift dependence of the characteristic frequency shifts is depicted. Solid line is 
i/*, and dashed line is z/j* pp . 
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We also simulate the line transfer with the cosmological parameters found by Spergel et 
al. (2003) in their analysis of WMAP data, who suggest h = 0.7ll°;^, Q M h 2 = 0.135±g;gg|, 
Q h h 2 Q = 0.0224±g;gg|, and Q A = 1.0±g;g| - ft M - Fig.8 presents the results at z s = 10. The 
solid line adopts the cosmological parameters of LR99, while the dashed line uses those of 
Spergel et al. (2003). We find the effect of the cosmological parameters on the line profile is 
small. 

1.2e-04 r- 
1.0e-04 - 
8.0e-05 - 

"i 

^ 6.0e-05 - 
4.0e-05 - 
2.0e-05 - 



1.33 1.34 1.35 1.36 1.37 1.38 1.39 1.4 

/L, [jU m] 

Fig. 8. — Observational line profiles at z s = 10 are presented. Solid line is the model with 
LR99 cosmological parameters, while dotted line is those of WMAP. 




5. SUMMARY 

In this work we extended LR99 by improving the numerical treatments of the Lya opac- 
ity and scattering processes. As part of these improvements, the effect of asphericity in the 
expanding IGM on the Lja profiles was examined. The results confirm LR99 very well when 
the aspherical deceleration is small. The effect of aspherical expansion on the line profile is 
reflected in a broadening of the FWHM and a shift of the peak intensity to longer wave- 
lengths over that in the spherical model. The effect of the redshift of reionization is reflected 
in the normalization of the frequency shift and escape distance from the source rather than 
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in the shape of the normalized line ratio. Our results were insensitive to small changes in the 
cosmological parameters, with little difference in the results found using the cosmological 
model of LR99 and that derived recently from WMAP data. To detect both properties will 
be important to find the initial stage of the formation of LSS. Fortunately, it requires a 
spectroscopic resolution of R ~ 100 — 1000, which is within the design characteristics for 
future observational projects like JWST. 

We are grateful to the referee for his/her advisable comment and excellent refereeing, 
which improve our paper very much in both content and presentation. We thank Profs. S. 
Inagaki, S. Mineshige, T. Totani for their continuous encouragement. 

REFERENCES 

Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 
Barkana, R. & Loeb, A. 1999, ApJ, 523, 54 
Breit, G., & Teller, E. 1940, ApJ, 91, 215 

Dey, A., Spinrad, H., Stern, D., Graham, J. R., & Chaffee, F. H. 1998, ApJ, 498, L93 

Gnedin, N. Y. & Ostriker, J. P., 1997, ApJ, 486, 581 

Gunn, J. E., & Peterson, B. A. 1965, ApJ, 142, 1633 

Hu, E. M., Cowie, L. L., & McMahon, R. G. 1998, ApJ, 502, L99 

Kamaya, H., Silk, J., 2002, MNRAS 332, 251 

Kogut, A., et al. 2003, ApJS, 148, 161 

Lin, C. C, Mestel, L., & Shu, F. H., 1965, ApJ, 142, 1431 

Loeb, A., & Rybicki, G. B. 1999, ApJ, 524, 527 (LR99) 

Miralda-Escude, J., Haehnelt, M. & Rees, M. J. 2000, ApJ, 530, 1 

Shchekinov, Iu., 1991, Ap&SS, 175, 57 

Spergel, D. N., et al, 2003, ApJS, 148, 175 

Spinrad, FL, Stern, D., Bunker, A., Dey, A., Lanzetta, K., Yahil, A., Pascarelle, S., & 
Fernandez-Soto, A. 1998, AJ, 116, 2617 



-18- 



Spitzer, L. J., k Greenstein, J. L., 1951, ApJ, 114, 407 

Weymann, R. J., Stern, D., Bunker, A., Spinrad, H., Chaffee, F. H., Thompson, R. L, & 
Storrie-Lombardi, L. J. 1998, ApJ, 505, L95 

Witt, A. N. 1977, ApJS, 35, 1 



This preprint was prepared with the AAS IATgX macros v5.0. 



